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ABSTRACT 

We present the results of deep observations of the Ophiuchus cluster of galaxies with INTEGRAL in the 3-80 keV band. We analyse 3 
Ms of INTEGRAL data on the Ophiuchus cluster with the IBIS/ISGRI hard X-ray imager and the JEM-X X-ray monitor. In the X-ray 
band using JEM-X, we show that the source is extended, and that the morphology is compatible with the results found by previous 
missions. Above 20 keV, we show that the size of the source is slightly larger than the PSF of the instrument, and is consistent with 
the soft X-ray morphology found with JEM-X and ASCA. Thanks to the constraints on the temperature provided by JEM-X, we show 
that the spectrum of the cluster is not well fitted by a single-temperature thermal Bremsstrahlung model, and that another spectral 
component is needed to explain the high energy data. We detect the high energy tail with a higher detection significance (6 Act) than 
the BeppoSAX claim (2a). Because of the imaging capabilities of JEM-X and ISGRI, we are able to exclude the possibility that the 
excess emission comes from very hot regions or absorbed AGN, which proves that the excess emission is indeed of non-thermal 
origin. Using the available radio data together with the non-thermal hard X-ray flux, we estimate a magnetic field B ~ 0.1 - 0.2 pG. 
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1. Introduction 

Clusters of galaxies are the biggest bound structures of the uni- 
verse, and, according to the hierarchical scenario of structure for- 
mation, the latest ones to form. They are filled by a hot ( 1 7 - 1 8 
K) plasma, called intra-cluster medium (ICM), and thus radiate 
in soft X-ray bands through thermal Bremsstrahlung. 

The most massive clusters should form by merging of 
smaller clusters, that w ould create sho ck waves when the ICM 
of the clusters merges dSarazi n 1999). Such events are prob- 
ably the major source of heating in the ICM. Models predict 
that a large population of relativistic electrons should be cre- 
ated during the merging event, which are expected to radiate 
through synchrotron emission in the radio domain, and through 
inverse-Compton (IC) scattering with th e cosmic microwave 
backg round (CMB) in the hard X-ray band (Ens slin & Bier mann 
1998). Another possible model involves a population of multi- 
TeV electrons that would radiate in hard X-rays through syn- 
chrotron emission dTimokhin et al.ll2004l) . While the extended 
radio emission fro m several clusters has been known for a 
long time (see e.g. iFeretti & Giovanninll d2007l) ). there is still 
no firm detection of the hard X-ray emission. BeppoSAX ob- 
servations report on the dete ction of a hard tail in the X -ray 
spectrum of at least Coma dFusco-Femiano et al.l |2004|) and 
Abell 2256 dFusco-Femiano etaL 20051). but these detection s 
are rather weak and controversial (Rossetti & Molendil 12004). 
Hence, confirmation of these results by other instruments is im- 
portant. 

The Ophiuchus cluster is the second brightest cluster in th e 
2-10 keV band. It is a nearby (z=0.028. |johnston et al.l dl98ll) ) 
rich cluster located in the direction of the Galactic Center (I - 
0.5°, b = 9.4°), with a very high plasma temperature (kT ~ 10 



keV). Using ASCA data. I Watanabe et alj d2001l) showed that the 
cluster is not dynamically relaxed and exhibits some regions 
with very high temperature (kT > 13 keV). Because these char- 
acteristics are very similar to the Coma cluster, they concluded 
that this cluster has also experienced a major merging even t 
in the recent past. In the radio domain, Johnst on" et al.l dl981) 
claimed that the source is associated with the steep-spectrum ra- 
dio source MSH 17-203. The identification of this radio source 
as a radio halo from the Ophiuchus cluster implies the presence 
of relativistic electrons, and hence predicts the presence of a non- 
therma l high-energy tail in the X-ray spectrum. lNevalainen et al.l 
(2004) searched for such an excess in a sample of nearby clus- 
ters with BeppoSAXfPDS. For the Ophiuchus cluster, their anal- 
ysis reveals a mean temperature of 9.1 + 0.6 keV, and a 2<x ex- 
cess at hard X-rays. However, at such a low significance level, 
the excess might as well be due to statistical fluctuations or in- 
strument systematic effects. Moreover, the PDS instrument on 
board BeppoSAX was non-imaging, so it was impossible to ex- 
tract any information on the morphology of the hard X-ray emis- 
sion. Therefore, the spectrum might be contaminated by ab- 
sorbed point sources. 
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In this paper, we present the results of deep (3 Ms) observa- 
tions of the Ophiuchus cluster with the IBIS/ISGRI and JEM-X 
instruments on board INTEGRAL, in the aim of investigating the 
presence of a high-energy tail. We present the specific method 
we have used to analyse ISGRI data. Finally, we use radio data 
to consider possible models for the non-thermal emission of this 
object. 
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2. Data Analysis 

2.1. ISGRI data analysis 

The I BIS/ISGRI instrument on board INTEGRAL dLebrun et alj 
l2003h is a wide-field (29° x 29°) coded-mask instrument sen- 
sitive in the 15-400 keV band. Its angular resolution (12 ar- 
cmin FWHM nominal) is of the same order of magnitude as the 
size of the core of the cluster (diameter ~ 10', IWatanabe et al.l 
J2001R Because of the large field-of-view (FOV) and the fact 
that INTEGRAL spends a significant part of its observing time 
in the Galactic Bulge region, the amount of data accumulated by 
the instrument on the Ophiuchus cluster is quite large, which al- 
lows us to reach a high signal-to-noise ratio for the cluster in the 
17-60 keV band. Our analysis covers 1580 pointings (or Science 
Windows, ScWs), including 1493 ScWs from public data and 87 
from the INTEGRAL AO-4 Key Programme, for a total observ- 
ing time of 3 Ms. However, the location of the source (9.3° from 
the Galactic Center and 14.5° from the very bright X-ray binary 
Sco X-l) makes difficult the imaging analysis of this region with 
ISGRI. Furthermore, the periodic shape of the IBIS mask causes 
the images gene rated by the standard O ffline Scientific Analysis 
software (OSA, ICourvoisier et alj d2003l) ) to present ghosts of 
bright sources which are not completely removed at some spe- 
cific positions in the deconvolved sky images. This adds impor- 
tant systematic errors to long exposure mosaics in the vicinity of 
bright sources. 

In a coded-mask instrument, the sky images are produced 
by the deconvolution of the shadow patterns cast by the sources 
in the FOV on the plane of the sky. When observing a single 
source, a large fraction of the detector (> 45%) is not illumi- 
nated by the source, which allows us to measure the level of the 
background in each pointing. The detector image, or "shadow- 
gram", produced when observing a field containing n sources 
is a superposition of the shadow patterns of all the individual 
sources, 

// 

s (*,y) = • PIF * x >y) + bB< < x ^ (!) 

where fi and PIFjix, y) are respectively the flux of the /th source 
and the shadow pattern (called "pixel illumination fraction") 
it casts, and B(x,y) is a "background map", i.e. a model of 
the background, in the pixel with coordinates (x,y). Therefore, 
bB(x,y) gives the total background model in the pixel (x,y). 
If the background was flat over the whole detector, the de- 
convolution process would not be sensitive to the background. 
Therefore, background maps take the deviations from the flat 
background into account. The properties of the coded mask 
method allow us to compute a background map by averaging 
the count rates of a pixel for many pointings, excluding the pix- 
els illuminated by a strong source. To take the variability of the 
background into account, we computed a specific background 
map for each revolution of INTEGRAL. 

To the present day, the exact shape of the IBIS mask is not 
completely understood, so the PIFs created by the standard OSA 
software are not perfectly correct. The occurence of even a hand- 
ful of wrongly modelled pixels can very strongly affect image 
reconstruction. The identification of these pixels is made diffi- 
cult by the fact that bright sources do also affect counts in ex- 
posed pixels, which should not be removed. In order to identify 
badly modelled pixels, we procede in the following way: we fit 
the model described in Eq. (Q~|i to the shadowgram and extract 
the background level and the flux of the sources present in the 



FOV. Then we remove bright sources and background from the 
shadowgram: 

n 

S'(x,y) = S (x, y)~Yjfi- PIF i(x> J) ~ bB{x, y), (2) 

i=l 

with error map given by error propagation. In the case of a per- 
fect PIF, the pixel distribution of the transformed detector image 
would be Gaussian, and hence 99.7% of the pixels would have a 
value within [fj. - 3cr,/i + 3cr]. Since the ISGRI detector contains 
~ 15'000 pixels, statistically we expect ~ 45 of them to show a 
deviation from the mean value larger than 3cr. However, we no- 
tice that ~ 150 pixels deviate from the mean value by more than 
3cr, so we conclude that a large fraction of these pixels are incor- 
rectly modelled. Then we perform a standard OSA 7.0 analysis 
ignoring these pixels. We also neglect the outer parts of the im- 
age (with off-axis angle > 14°), where the systematic effects are 
the most important. Finally, we make a mosaic from the modi- 
fied images in a standard way. 

Unlike the image reconstruction process, the standard OSA 
spectral extraction procedure does not suffer from these prob- 
lems as long as the input source catalog contains all sources 
brighter or comparable with the analysed source. Therefore, for 
spectral extraction we performed a standard OSA 7.0 analysis 
with a source catalog extracted from the mosaic image. 

2.2. JEM-X data analysis 

The J EM-X X-ray monitor on board INTEGRAL dLund et al.l 
2003) consists of two identical X-ray detectors with coded mask, 
JEM-X 1 and JEM-X 2, sensitive in the 3-35 keV band. It is de- 
signed for the spectroscopic and imaging study of the sources 
detected by IBIS, with a better spatial resolution (3.35 arcmin 
FWHM nominal) and a field-of-view of 7.5° half-response in di- 
ameter. On the Ophiuchus cluster, we reached a 160 ks effective 
exposure time, for a total observing time of 350 ks. 

The analysis of the Ophiuchus field with JEM-X does not 
present the same difficulties as those affecting the IBIS obser- 
vations. Indeed, both the Galactic Center and Sco X-l are out- 
side the field of view of JEM-X. We therefore performed a com- 
pletely standard image analysis, taking advantage of the im- 
provements in the image reconstruction process introduced in 
OSA 7.0, in particular regarding the astrometry. In order to reach 
the deepest sensitivity of JEM-X, we combined all JEM-X 1 and 
JEM-X 2 ScWs in a single mosaic. Spectral extraction has been 
performed by extracting the fluxes from this mosaic in four sep- 
arate energy bands chosen to provide clear detection in each in- 
dividual band and to complement the ISGRI spectrum without 
gap. The fact that the Ophiuchus cluster is probably extended 
with JEM-X is not a problem, since the flux is integrated over the 
source. Therefore we simply let the width of the fitted Gaussian 
free. 

When extracting fluxes from JEM-X mosaics, there is no 
standard way in OSA to provide associated responses. The re- 
distribution matrices (RMFs) of both JEM-X instruments are sta- 
ble with time and quite close for both instruments; we therefore 
simply use the RMF of JEM-X 1. The area responses (ARFs) 
are however varying both with time and with the instrument. In 
particular, the low-energy efficiency changed very strongly (by 
a factor close to 3 for JEM-X 1) when the event rejection cri- 
teria were modified. It is therefore necessary to build an ARF 
suitable for the analysis of the spectrum extracted from the com- 
bined mosaic. To do this, we extract the correct ARF for each 
ScW and calculate the average of these ARFs weighted by the 
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Fig. 1. JEM-X significance image in the 3-18 keV band, with 
surface brightness contours from ASCA overlayed. The inset 
shows the image of a known point source in the same field, V 
2400 Oph, to show the extended nature of the source. The cross 
shows the position of 1RXS J171209.5-231005, the nearest X- 
ray point source. 

effective exposure time for each ScW, taking into account the 
dead time and the vignetting. As most observations use JEM-X 
1 with the new rejection criterion, in practice the resulting ARF 
is quite close to that valid for any of the late JEM-X 1 ScWs. 
Therefore, the choice of a particular ARF does not affect quali- 
tatively any of the results presented in this paper. 

3. Imaging results 

As demonstrated in Eck ert etaLI d2007) for the Coma cluster, the 
angular resolution of INTEGRAL can be sufficient to get impor- 
tant information on the morphology of the hard X-ray emission 
from clusters. Since the Ophiuchus cluster is a very bright clus- 
ter below 10 keV, one can use the JEM-X monitor to study the 
change in morphology with energy. Fig.Q]shows the JEM-X mo- 
saic image of the cluster in the 3-18 keV band, with contours 
from ASCA/GIS observations of the cluster. For comparison, the 
inset shows the JEM-X image of the cataclysmic variable V 2400 
Oph, in the same field. The extended nature of the Ophiuchus 
cluster in the JEM-X images is clear. 

Because we detect only the core of the cluster, it is not pos- 
sible to fit the surface brightness with an isothermal beta profile, 
but fitting the image with a Gaussian model, 

/(r)=Aexp|-ln(2)^V (3) 

where R is the half-width at half-maximum (HWHM), one can 
get an indication of the size of the detected region and of possi- 
ble deviations from the spherical model. Fitting the JEM-X im- 
age shown in Fig. Q] with the model described in Eq. (01, we 
find R = 3.6 + 0.1 arcmin. This value is significantly larger 
than the PSF of JEM-X (1.8 arcmin HWHM), which confirms 
that the source is s een as extended. This is, after SN 1006 
dKalemci et al. 1 120061) . the second report of an extended source 
for JEM-X. The R parameter corresponds to the superposition of 
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Fig. 2. Radial profile of the source with a core radius of 3.1' 
extracted from JEM-X data (dashed blue line), compared to the 
PSF of the instrument (dashed green). The solid red line shows 
the total apparent JEM-X profile. The purple bars show the data 
points extracted from the JEM-X image. 



the radial profile of the source and of the PSF of the instrument, 
^ 2 = ^source + H Mjj^ x . Therefore, the angular size of the 
source at half-maximum becomes R SOU rce = 3.1 +0.1 arcmin. 
Figure [2] shows the radial profile of the source compared to the 
PSF of the instrument, as well as the total apparent radial pro- 
file from JEM-X data. We can see that the radial profile is well 
modelled by a Gauss ian. The extensio n of the source is clear. 

From ASCA data, Watana be et al.l(l2001l) find the same value 
for the core of the cluster, R c = 3.2'. They also report on an ex- 
tended region North- West of the center of the cluster that shows 
significant deviations from the spherical model. The authors note 
that this excess coincides with the closest X-ray point source, 
1RXS J171209.5-231005. Figure shows the residuals of the 
JEM-X image compared to the spherically symmetric model. A 
clear excess North- West of the center is also visible. However, it 
is clear that this excess is not due to 1RXS J171209. 5-231005, 
which is located outside the region detected by JEM-X. The de- 
viations from the spherical model are hence more likely due to 
over-densities or excess temperatures in the ICM, which indi- 
cates that the cluster is not dynamically relaxed, probably be- 
cause of a merger in its recent history. 

We used the available JEM-X data to create a hardness ratio 
map (7-18 keV to 3-7 keV) in order to estimate the temperature 
variations within the detected region of the cluster. To estimate 
the dependence of the JEM-X 7-18 to 3-7 keV hardness ratio 
with temperature, we simulated JEM-X spectra for models with 
different input temperatures and computed the hardness ratios. 
The result is shown in the top panel of Fig. [4] Because of the 
sensitivity up to 18 keV, the temperature dependence is steeper 
than in the case of ASCA. In the center of the cluster, we find a 
hardness ratio HR = 0.78 + 0.09, which corresponds to a tem- 
perature kT = 9.1+M. The bottom panel of Fig. [4] shows the 
deviations of the hardness ratio compared to the central value, in 
units of a. We find no deviations to the central value at signif- 
icance level above 0.7<x, which corresponds to variations of +2 
keV compared to the central value. In ASCA data, although there 
is evidence for a strong spatial dependence of the temperature, 
this dependence is found only in the outer parts of the cluster, 
while the core, which is the only region detected by JEM-X, 
shows an almost iso-thermal profile, compatible with the JEM- 
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Fig. 3. Residuals from the spherically symmetric model with 
JEM-X brightness contours overlaid, in units of <x. Similar to 
previous ASCA observations, a significant excess is found North- 
West of the center of the cluster. 

X results. This analysis shows that a single-temperature model 
should describe well the total JEM-X spectrum. 

Fig. [5]shows the ISGRI significance image in the 20-40 keV 
band extracted from the mosaic created with the method de- 
scribed in Sect. I2.ll Because of the lower spatial resolution, the 
extension of the source is not obvious (see the image of V 2400 
Oph in the inset for comparison). However, a more detailed anal- 
ysis of the source reveals that the source is indeed extended also 
for ISGRI. Fitting the image with the model described in Eq. 
0), one finds R = 8.9 + 0.3 arcmin, which is slightly larger 
than the PSF of the instrument (7.2 arcmin HWHM). Taking into 
account the contribution of the radial profile of the source and 
of the PSF of the instrument similar to the analysis presented 
above for JEM-X, this corresponds to R S0U rce = 5.1*g5. Even 
though the source is seen as extended, the apparent size of the 
Ophiuchus cluster is smaller than the PSF of the instrument. To 
estimate the error on the standa r d flux extraction in the case of 
extended sources, iRenaud et al.l ((2006) presented the results of 
simulations assuming a uniform disk. For an angular size of 5', 
they evaluate the relative error to be ~ 5%. Since the cluster has 
a peaked radial profile instead of a uniform disk, the relative er- 
ror on the flux extracted by the OSA spectral extraction tool is 
< 5%. For this reason, from now on we will treat the source as 
point-like for ISGRI. 

Since the apparent size of the source is smaller than the PSF 
of ISGRI, it is not possible to study in detail the hard X-ray mor- 
phology of the source. However, it is possible to measure the 
best fit position for the center of the source in different energy 
bands. Figure [6] shows the 90% error circles for the position of 
the source in the 20-24, 24-30 and 30-40 keV bands overlaid on 
the JEM-X 3-18 keV mosaic image. The yellow cross shows the 
position of the center of the cluster from the Chandra image. 
The black crosses represent the position of the 2 brightest point 
sources close to the center in the Chandra image. We can see 




0.55 1 1 1 1 1 1 1 1 

6 7 8 9 10 11 12 13 
kT [keV] 




2 -1.5 -1 -0.5 0.5 1 1.5 



Fig. 4. Top: Simulated JEM-X 7-18/3-7 keV hardness ratio as 
a function of temperature. Bottom: Deviations of the JEM-X 7- 
18/3-7 keV hardness ratio compared to the central value in units 
of cr, with JEM-X brightness contours overlaid. The map does 
not show any significant spatial dependence of the temperature 
within the detected region. 



10' 




5 10 15 20 



Fig. 5. ISGRI significance map in the 20-40 keV band, with sur- 
face brightness contours from ASCA in green. The inset shows 
the ISGRI image of a point source, the cataclysmic variable V 
2400 Oph, for comparison. 
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on this image that the hard X-ray emission is compatible with 
the soft X-ray morphology of the cluster, and is not displaced to- 
wards any of the weak point sources seen in the Chandra image, 
which excludes the possibility that these point sources contribute 
in a significant way to the high-energy spectrum. 
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Fig. 6. 90% error circles for the position of the center of the 
source in the 20-24 (red), 24-30 (blue) and 30-40 keV (white) 
ISGRI mosaic images, overlaid on the JEM-X 3-18 keV image. 
The yellow cross shows the position of the centroid of the cluster 
from the Chandra image. The black crosses show the position of 
the two brightest X-ray point sources detected in the Chandra 
image. 



4. Spectral analysis 

4. 1 . INTEGRAL broad band spectrum of the cluster 

Since, unlike the case of the Coma cluster, the source is only very 
slightly larger than the PSF of ISGRI, we can use the standard 
OSA 7.0 spectral extraction tool with our time-dependent back- 
ground maps to get an accurate flux estimation. Starting from 
OSA version 7.0, the calibration of ISGRI is now valid down 
to 17 keV. Because of the large exposure time, we were able 
to reach a high signal-to-noise ratio in both JEM-X (3-18 keV 
band) and ISGRI (17-60 keV band). In the 40-60 keV band, the 
detection significance reaches 6cr. In the 60-80 keV band, the 
source is marginally detected by ISGRI at the 1.5<x level. It is 
important to note that no cross-calibration factor between ISGRI 
and JEM-X data was needed for this work. Introducing such a 
factor in the XSPEC fitting model results in a cross-calibration 
of 1 .0, which shows the excellent calibration work done by the 
ISGRI and JEM-X teams for the release of OSA 7.0. 

Fitting the combined JE M-X/ISGRI data with a single- 
temperature MEKAL model dKaastra & Mewdl2000l). with the 
abund ance fixed to 0.49 compared to t he solar val ue ( Moh r"et~aT] 
119991) and redshift fixed to z = 0.028 dJohnston et all 198 lb . one 
finds a bad representation of the data, with a reduced^ 2 exceed- 
ing 2.3. Another spectral component is hence needed to get a 
better fit to the data. Fig. [7] shows the combined JEM-X/ISGRI 



Table 1. Results of the fitting procedure for a thermal + power- 
law model, with abundance and redshift values fixed to the litter- 
at ure value and phot o n inde x fixed to a - 2.0, similar to the work 
of iNevalainen et all d2004l) . In Method 1, the temperature and 
the normalisation of the Bremsstrahlung component are left free 
while fitting. In Method 2, we used the data in the 3-20 keV band 
to fix the temperature and normalisation of the thermal compo- 
nent, and let free only the normalization of the non-thermal com- 
ponent. a PF-test is the null hypothesis probability when adding 
the non-thermal component given by the F-test. fe 20-60 keV flux 
of the non-thermal component, in units of 10" 12 ergs s crrT 2 . 
The errors are quoted at the lcr level. c Confidence level for the 
detection of the non-thermal component. 
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Fig. 7. JEM-X/ISGRI combined spectrum in the 3-80 keV band. 
The solid line is a fit to the 3-20 keV part of the spectrum with 
a single-temperature MEKAL model. The bottom panel shows 
the residuals from the model. The hard X-ray excess is clear. 



spectrum and the residuals from the thermal model with tem- 
perature kT = 8.50 keV fitted to the spectrum below 20 keV. 
The model exhibits significant residuals above 20 keV, which 
could be due to th e presence of an addition al non-thermal com- 
ponent. Similar to Neva lainen et al.l d2004). we assume that the 
non-thermal component is a power-law with a fixed photon index 
a = 2.0. Fig. |8] shows the unfolded spectrum with the thermal 
model and an additional power-law component fitted simultane- 
ously (Method 1). Table[T]shows the results of the fitting proce- 
dure using such a power-law in addition to the thermal emission, 
with the plasma temperature let free while fitting (Method 1) and 
fixed to the value found by fitting independantly the 3-20 keV 
part of the spectrum (Method 2), which is expected to be com- 
pletely dominated by the thermal component. The BeppoSAX 
claim (2cr excess at high energies) was obtained by fixing the 
temperature to the value obtained from the low-energy spectrum, 
which is similar to Method 2. Therefore, we can see that the clear 
excess (6.4<x) detected with this method in INTEGRAL data is 
much more significant than the value obtained by BeppoSAX. 
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Fig. 8. INTEGRAL unfolded spectrum of the Ophiuchus cluster 
in the 3-80 keV band fitted by a MEKAL+power-law model with 
fixed photon index a = 2.0 (green solid line). The red and green 
dashed lines represent the contribution of the thermal, respec- 
tively non-thermal components as a function of energy. 



The value found for the plasma temperature is consistent 
with the BeppoSAX measurement of kT = 9.1 + 0.6 keV, but 
significantly lower than the ASCA value of kT = 10.9 + 0.4 keV. 
However, the ASCA/GIS instrument was sensitive only up to 10 
keV, and thus could not detect the cutoff. Because of the sensitiv- 
ity above 10 keV, the INTEGRAL and BeppoSAX values certainly 
provide more accurate measurements. 

Fixing the thermal component to the value found in the 3- 
20 keV band, one can try to constrain the acceptable range for 
the photon index of the non-thermal component. Leaving only 
the photon index and the normalisation of the power-law free 
while fitting, one finds a - 1.62*?g?. The photon index of 2.0 
used to derive the values of the HXR emission is compatible 
(within lcr) with this value. With such a hard photon index, it 
is clear that the excess cannot be explained by the presence of 
hot gas. Indeed, if instead of a power-law we fit the excess with 
a second thermal component, we can derive a lower limit to the 
temperature of the second component, which we find to be kT > 
50 keV. Together with the temperature map from both JEM-X 
and ASCA, this proves that the excess cannot be explained by 
the presence of very hot gas. 



4.2. Possible origins of the non-thermal component 

Since the hard X-ray flux observed by INTEGRAL cannot be ex- 
plained by thermal Bremsstrahlung, there are two possible ori- 
gins for the non-thermal flux in the hard X-ray band. Indeed, 
it can be due either to in yerse-Compton (IC) scatter ing of GeV 
electrons with the CMB (Ensslin & Biermann 1998), or to syn- 
chrotron emission from an other population of very high energy 
electro ns (E ~ PeV, e.g. iTimokhin et al.1 (120041) . Ilnoue et all 
(2005;)). To construct our models, we assume that the electron 
distribution is a power-law with spectral index 2.0 and a high- 
energy cut-off. 

Assuming that the non-thermal component comes from IC 
scattering of the same electrons that produce radio halos, one 
can try to constrain the phys ical parameters nee ded to explain 
the radio/hard X-ray flux dThierbach et al.ll2003l) . in particular 
the value of the magnetic field in the cluster. Unfortunately, the 
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Fig. 9. Models for the spectral energy distribution of the non- 
thermal emission of the cluster, for a population of E ~ PeV 
electrons (top) and of E ~ GeV electrons (bottom), for differ- 
ent values of magnetic field and cut-off energy. The solid black 
line shows the 10 mCrab flux level in the TeV range. The radio 
data are from Joh nston et al.l d 198 lb . The dashed lines at ener- 
gies ~ 10 TeV show the attenuation of the IC spectrum due to 
absorption by the extragalactic infrared and CMB backgrounds. 



existing radio data are old, and the size of the region used to ex- 
tract the radio flux differs for all the measurements, so the shape 
of the radio spect rum cannot be well co nstrained. Using the radio 
data presented in Johnst on" et al.l (Il98ll) . we tried nevertheless to 
consider jointly the radio and hard X-ray data of the cluster to 
constrain models for the spectral energy distribution. The bot- 
tom panel of Fig.|9]shows the corresponding spectral energy dis- 
tribution, with models computed for several values of magnetic 
field B and cut-off energy E cut . 

Because of the low quality of the radio data, we cannot con- 
strain the properties of the electron population. However, the ra- 
tio between the synchrotron and IC components is given by 



R sync Mb 



IC 



B 2 /8n 

"CMB 



(4) 



The only unknown factor in this ratio is the magnetic field, so 
we can use the mean value of the radio data to estimate the value 
of B. Figure|9]shows the spectral models computed for different 
values of B. We can see that the best estimate for the magnetic 
field is 



B ~ 0.1 -0.2 /zG. 

This value is similar to 
X-ray radiation of the 



the value obtained 
Coma cluster (B 



from 



(5) 

the hard 
0.15 pG, 

Fusco-Femiano et al. (1999)), but contradicts the much higher 
values obtained from Faraday rotation measur es in a large sam- 
ple of radio-emitting clusters dKim et al ll99lh . 

If the hard X-ray emission comes from synchrotron radiation 
from a different population of electrons with highly relativistic 
velocities, we expect that the peak of the IC emission will be in 
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the TeV domain. The top panel of Fig.[9]shows the correspond- 
ing spectral energy distributions for different values of magnetic 
field and a cut-off energy E cut = 0.3 PeV. The black line shows 
the 10 mCrab flux level in the TeV domain. We can see that the 
magnetic-field values deduced from Faraday rotation measures 
(~ 1 — 5 pG) predict a TeV flux which is firmly detectable by the 
present generation of Cherenkov telescopes. 

5. Discussion 

The combined ISGRI/JEM-X spectrum of the cluster exhibits a 
clear {6 Act) excess at high energies. At such significance level, 
it is now clear that an additional component is needed to explain 
the relatively high hard X-ray flux observed both by INTEGRAL 
and BeppoSAX. Unlike the PDS instrument on board BeppoSAX, 
the ISGRI and JEM-X instruments on board INTEGRAL are 
imaging, which allowed us to study the changes in temperature 
throughout the detected region, and show that the hard excess 
cannot be explained by the presence of very hot thermal plasma 
in some regions of the cluster. Besides, in the large FOV of the 
PDS instrument, there was no guarantee that the possible excess 
was not the result of isolated point sources, for instance absorbed 
active galactic nuclei. The imaging capabilities of INTEGRAL 
allow us to exclude this possibility. Indeed, we find evidence 
that the source is extended even above 20 keV, and coincides 
spatially with the X-ray emission. Moreover, high angular reso- 
lution X-ray imaging of the cluster with Chandra does not reveal 
the presence of any bright point source below 10 keV in the re- 
gion detected by INTEGRAL, which shows that contamination 
of the INTEGRAL spectrum by point sources is very unlikely. 
Therefore, we claim with good confidence that the high-energy 
tail observed by INTEGRAL and BeppoSAX is certainly due to 
non-thermal emission from relativistic particles accelerated dur- 
ing a major merging event in the recent past. 

Assuming that the non-thermal component is explained by 
IC scattering from the same electrons that produce the syn- 
chrotron radio halo, we used the existing radio data together with 
the INTEGRAL measurement to make a joint model for the syn- 
chrotron/IC components, and estimated the magnetic field value 
to be B ~ 0. 1 -0.2 pG. This value is similar to that estimated with 
the sa me method for the Coma cluster (iFusco-Femiano et al.l 
1 1999b and with the typic al magnetic field value s extracted from 
equipartition assumption (Thierbach et al. 2003), but contradicts 
the higher valu es (B ~ 1 - 5 pG) derived from Far aday rota- 
tion measures dKim et al.l d 19911) . IClarke et ail (|2001|) ). To re- 
solve the discrepanc y between the Faraday ro t ation a nd Inverse- 
Compton methods, Goldshmidt & Rephaelj] d 1993b suggested 
that the magnetic field might be decreasing with radius, such 
that in the outer regions the electron cooling through IC emis- 
sion would dominate, while in the center, the electrons would 
radiate more strongly in the radio domain. However, the spatial 
resolution of INTEGRAL allows us to measure the magnetic field 
within a region of radius less than 6 arcmin from the center of 
the cluster, and therefore to show that the value of B measured 
from the HXR radiation does not come from th e outer regions 
of the cluster. More recently, iBeck et all (|2003) suggested that 
the presence of turbulence might introduce a bias in the value of 
B derived from Farad ay rotation measures. F rom X-ray imaging 
of the Coma cluster, Schuec ker et al.l d2004t) demonstrated that 
turbulence is indeed playing an important role in unrelaxed clus- 
ters. Although it is much less studied, Ophiuchus is similar to 
Coma in many aspects, so this explanation could apply as well 
in our case. In any case, the measurement of the magnetic field 
through Faraday rotation measures relies on several assump- 



tions, whereas the expected radio/hard X-ray flux correlation 
provides a direct measurement of the magnetic field, so the lower 
values obtained in both Ophiuchus and Coma (B ~ 0.1 pG) cer- 
tainly provide a better estimate of the cluster magnetic field. 
Better constraints on the shape of the radio spectrum are there- 
fore important to measure the magnetic field with better accu- 
racy. 

On the other hand, if the magnetic field values derived from 
Faraday rotation measure are correct, the IC emission from the 
population of electrons that radiate in the radio domain is not 
sufficient to explain the hard X-ray flux detected by INTEGRAL. 
This implies the presence of another population of electrons at 
much higher energies (E ~ 100 TeV) that would radiate in hard 
X-rays through synchrotron radiation and would up-scatter the 
photons of the CMB to produce TeV emission (see the top panel 
of Fig. |9). In this case, our models predict a flux in the TeV do- 
main that would be firmly detectable by HESS. Indeed, our low- 
est curve, which corresponds to the highest allowed value for 
the magnetic field (10 pG), coincides with the sensitivity limit of 
HESS. Therefore, observations of the cluster in very high ener- 
gies are crucial to constrain the models of particle acceleration 
in Ophiuchus in particular and clusters of galaxies in general. 

6. Conclusion 

In this work, we presented the results of broad-band INTEGRAL 
observations of the Ophiuchus cluster in the 3-18 keV band with 
the JEM-X X-ray monitor and in the 17-80 keV with the ISGRI 
hard X-ray imager. We have shown that the source is extended 
for both instruments and that the morphology of the source is 
consistent with previous studies (Sect.O. We presented a JEM- 
X hardness-ratio map (7-18 keV to 3-7 keV) and found that the 
plasma temperature does not vary significantly within the de- 
tected region. In Sect. [4] we presented the total INTEGRAL spec- 
trum of the cluster in the 3-80 keV band. We explained that the 
emission cannot be described by a single-temperature MEKAL 
model, and that an additional spectral component is needed to 
explain the significant excess found at energies above 20 keV 
compared to the best fit (kT = 8.5 keV) for the thermal emis- 
sion below 20 keV. Consequently, we measure a hard X-ray ex- 
cess emission from the Ophiuchus cluster, at a much higher con- 
fidence level (6.4cr for the fixed temperature model) than the 
BeppoSAX claim. 

With spectral and imaging analysis, we have shown that this 
excess cannot be explained by the presence of very hot ther- 
mal gas or by contamination of the spectrum by point sources. 
Indeed, fitting the excess with a second thermal component, we 
found an unrealistic lower limit to the temperature, kT > 50 
keV. Comparing the JEM-X and ISGRI images with X-ray im- 
ages from higher resolution instruments (Chandra and ASCA), 
we found that the morphology of the cluster does not depend 
on energy, and that no bright point sources are detected in the 
high-resolutions images, which excludes the possibility that the 
hard X-ray emission comes from point sources embedded in the 
cluster. Therefore, we conclude that the high-energy tail is due to 
non-thermal emission from relativistic electrons, probably accel- 
erated during a merging event in the recent history of the cluster. 

To investigate the origin of this non-thermal component, we 
presented different models for the spectral energy distribution 
(Sect. 14.2b . If the emission is of inverse-Compton origin from 
GeV electrons on the CMB, we used the available radio data to 
construct a joint radio/hard X-ray model for the spectral energy 
distribution. Comparing the level of the radio and hard X-ray 
flux, we were able to estimate the magnetic field in the cluster, 
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which we found to be B ~ 0.1 - 0.2 ^G. On the other hand, if 
the hard X-ray emission comes from synchrotron radiation of 
another population of very high energy electrons, our models 
predict a detectable flux in the TeV domain. Therefore, obser- 
vations of the Ophiuchus cluster in the TeV domain, as well as 
better measurements of the radio spectrum of the cluster, will 
allow us to constrain the models, and determine if clusters of 
galaxies are able to accelerate particles up to very high energies. 
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